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Abstract. We present recent advances on Dirichlet forms methods either to extend fi- 
nancial models beyond the usual stochastic calculus or to study stochastic models with 
Q_i ■ less classical tools. In this spirit, we interpret the asymptotic error on the solution of an 
sde due to the Euler scheme (Kurtz and Protter [Ku-Pr-91a]) in terms of a Dirichlet form 
on the Wiener space, what allows to propagate this error thanks to functional calculus. 
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t^. ! Introduction 

Considering a Dirichlet form amounts to consider a strongly continuous symmetric 
contraction semi-group on an L 2 -space which possesses in addition the property of being 
positive on positive functions (cf. [Fu-Os-Ta-94], [Bo-Hi-91 ], [Ma-R6-92]). It is a par- 
ticular case of Markovian potential theory, with several special features due to the use of 
Hilbertian techniques and to the fact that positivity and contraction properties extend to 
infinite dimensional framework thanks to Fatou's lemma in measure theory. Many Dirich- 
let structures are constructively obtained on the Wiener space and on the fundamental 
^ ■ spaces of probability theory (Poisson space, Monte Carlo space) which may be thought as 
hypotheses in order to study error propagation through stochastic models (cf. [Bou-03b]). 

Since the discovery by M. Fukushima, at the end of the seventies, that Dirichlet forms 
allow to extend the stochastic calculus to processes which are not semi-martingales (cf. 
[Fuk-80]) a lot of works have been developed in this direction, even beyond the Dirichlet 
forms framework. To this extend we quote the approach to time-dependent Dirichlet 
forms developed by Oshima [Osh-92] and the more recent approach of Stannat [Sta-99] 
and Trutnau [Tru-OO] about a new theory of generalized Dirichlet forms. As in finance the 
heart of the complete market property and more generally of the portfolio management 
is the stochastic integral, a particular interest has been devoted to methods giving rise to 
new stochastic integrals. 

We shall give, at first, a short outlook on recent results related to Dirichlet forms and 
connected with financial motivations. We include some Malliavin calculus approaches 
when they amount to the use of the Ornstein-Uhlenbeck structure on the Wiener space. 
After recalling, in a second part, the main properties of Dirichlet forms and the interpreta- 
tion of the functional calculus on the squared field operator in terms of error propagation, 
we focuse, in a third part, on the question of the asymptotic error due to the resolution 
of a stochastic differential equation by the Euler scheme. We show that the asymptotic 



to propagate the error on the example of a level volatility model for pricing and hedg- 
ing procedures. We put the general question of the validity of such a propagation as an 
asymptotic calculus principle, and we give partial arguments for this principle. 

I. Some recent works. 

First must be mentionned the idea of using Malliavin's integration by parts technique 
to speed up the computation of the Greeks or other quantities in finance. After the 
collective papers of Fournier and al [Fo-La-Le-Li-To-99] [Fo-La-Le-Li-01], improvements 
have been brought to complex options [Go-Ko-01] and to the more general question of the 
sensitivity to some parameters with the aim of calibration of a model. As integration by 
parts formulae are available in more general Dirichlet forms situations than the Ornstein- 
Uhlenbeck structure on the Wiener space (cf. [Bou-03b] Chapter V), the same approach 
may be performed for instance on the Poisson space for studying models with jumps 
[El-Pr-04]. 

One of the first success of Malliavin calculus was about proving existence of densities 
for solutions of sde's with smooth coefficients and Dirichlet forms methods have been able 
to extend such results to the case of Lipschitz coefficients [Bo-Hi-91]. Several authors 
remarked that these tools give also means of improving the computation of densities and 
establishing estimates for the laws of random variables with some regularity assumption. 
Let us quote ([Ko-Pe-02] , [Ca-Fe-Nu-98] , [Bo-Ek-To-04], [Bou-05b]) whose results aren't 
limited to applications in finance. With suitable hypotheses it is possible, to get explicit 
closed formulae for the density even with some liberty in the choice of a weight function 
allowing an optimization for Monte Carlo simulation. 

After the classical works of M. Fukushima and Y. Le Jan on stochastic calculus for ad- 
ditive functionals of symmetric Markov processes associated with a Dirichlet form [LeJ-78] 
the role of past and future a-algebras have been clarified by Lyons and Zheng (cf [Ly- 
Zh-98] [Tru-00]) and the main current of research, in order to leave the semi- martingale 
context, starts with the abstract definition of a Dirichlet process as sum of a local mar- 
tingale and a process with zero quadratic variation (see [F61-80]). Because the quadratic 
variation, as formal Dirichlet form, does not possess the closedness property, the Dirichlet 
form framework is replaced here by functional analytic arguments. The integral is gen- 
erally defined by a discretization procedure (cf. [F61-81], [Bou-85], [F6-Pr-Sh-95] ) or by 
a regularization procedure (see [Ru-Va-95], [Ru-Va-96]). These ways have been deepened 
with the center example of the fractional Brownian motion (cf. [Er-Ru-98], [Zah-98], 
[Fe-LaP-99], [Al-Ma-Nu-00] , [Ru-Va-00], [Gr-No-03], [Gr-Ru-Va-03]). The connections of 
these works with finance are many : attempting to generalize Girsanov theorem in or- 
der to define martingale measures by erasing more general drifts and using generalized 
stochastic integration (forward, symmetric and backward integrals) in order to deal with 
exotic models (cf [Fl-Ru-Wo-03]). About "inside trading" and the use of forward integral 
it is worth to quote [Le-Na-Nu-03]. 

At last, let us mention some uses of Dirichlet forms or Malliavin calculus to deal with 
processes with jumps by equipping the general Poisson space with a differential structure 
(cf. [De-Gr-Po-99], [Me-Pr-03], ) and the forthcoming book of P. Malliavin and A. Thal- 
maier [Ma-Th-05] whose last chapter is devoted to calculus of variations for markets with 
jumps, the other ones being strongly related with the above topics. 



II. Dirichlet forms theory seen as error propagation theory. 



we consider an erroneous quantity with a centered small error and apply to it a non linear 
map, we observe by an easy Taylor expansion argument that 

- the error is no more centered in general : a bias appears 

- the variance transmit with a first order calculus. 

Now if we go on, applying anew several non-linear applications 

- the variances and the biases keep (except special cases) the same order of magnitude 

- the biases follow a second order differential calculus involving the variances. 
With natural notation 

a n+l = fn+li x n) a n 

bias n+ i = /^ +1 (x„)bias„ + ^fn +l {x n )a 2 n . 

The first relation has been discovered, even in several dimension with correlation between 
the errors, by Gauss at the beginning of the nineteenth century. 

From this observation, in order to represent the propagation of small errors we may 
consider that 

1) the variances of errors have to be managed by a quadratic first order differential 
operator T, 

2) the biases of errors have to be represented by a linear second order differential 
operator A, 

the propagation of errors being the result of the following change of variable formulae : 
T[F(X 1 , . . . , X m ), G(Y U . . . , Y n )\ = ^ F'iXi, . . . , XJG'^Y,, F n )r[X 4 , YA 

ij 

A[F(X 1 , . . . , X m )\ = ]T F!(X U . . . , X m )A[X t ] + X - ]T F%(X U . . . , X m )r[X,, XA. 

i ij 

Because of these propagation rules for the variances and the biases, little errors may be 
thought as second order vectors. This old notion of differential geometry has been revived 
at the beginning of the eighties by the study of semi- martingales on manifolds (cf. [Sch-82] 
[Mey-82] [Eme-89]). 

Now, instead of germs of semi- martingales and second order vectors, we will use Dirich- 
let forms, carre du champ and generator. There are two important reasons for this, that 
I shall give just after recalling some definitions and examples. 

Definition An error structure is a term 

S= (n,A,T,D,T) 
where (Q, A, P) is a probability space, and: 

(1) ID is a dense subvector space of L 2 (Q, A,P) (also denoted L 2 (JP)); 

(2) r is a positive symmetric bilinear application from ID x ID into ^ 1 (1P) satisfying 
"the functional calculus of class C 1 fl Lip". This expression means 

\/u G E> m , W G E> n , VF: lR m -> H, VG: HT -> H 

with F, G being of class C 1 and Lipschitzian, we have F(u) G ID, G(v) G ID and 



(3) the bilinear form £[u,v] = ^JE[T[u,v}] is "closed". This means that the space ID 
equipped with the norm 

1 /2 

IMIb = (|M|£ 2(P) +£[u,ufj 

is complete. 
If, in addition 

(4) the constant function 1 belongs to ID (which implies T[l] = by property 2), we say 
that the error structure is Markovian. 

We will always write £[u] for £[u,u] and T[u] for F[u, u\. 

With this definition, the form £ is known in the literature as a local Dirichlet form on 
L 2 (fl, A, P) that possesses a "squared field" operator (or a "carre du champ" operator) 
T. These notions are usually studied on a- finite measurable spaces. We limit ourselves 
herein to probability spaces both for the sake of simplicity and because we will use images 
and products of error structures. 

Under very weak additional assumptions, to an error structure (also to a Dirichlet form 
on a (T-finite measurable space) a strongly-continuous contraction semigroup (P t ) t>0 on 
L 2 (P) can be uniquely associated, which is symmetric with respect to P and sub-Markov. 
This semigroup has a generator (A, DA), a self-adjoint operator that satisfies: 

clF 1 f) 2 F 

A[F(u)} = E ^-(^M + g E ^» r K u 3 \ P-a.s. 

% i i,j * 3 

for F: P m — > P of class C 2 with bounded derivatives and u G (VA) m such that r[tij] G 
L 2 (P). 

Example l.(Ornstein-Uhlenbeck structure in dimension 1) 

Q, = IR, ^ = Borel a-field £(P), P = W(0, 1) reduced normal law, D = 
/^(A/XO, 1)) = {« 6 L 2 (P),w' in the distribution sense belongs to L 2 (JP)}, T[u] = u' 2 , 
then (P, £>(P), A/"(0, 1), if 1 (A/'(0, 1)), T) is an error structure with generator 

PA = {/ G L 2 (P): /" - xf in the distribution sense G L 2 (P)} 

4f = -/" --/•/' 
2 2 

where / is the identity map on P. 

Example 2. (Monte Carlo structure in dimension 1) 

f2 = [0, 1], A = Borel a-field, P = Lebesgue measure, D = {hG L 2 ([0, l],dx) 
the derivative u' in the distribution sense over ]0, 1[ belongs to L 2 ([0, 1], c&r)}, r[iz] = u' 2 . 
Example 3. (Friedrich extension of a symmetric operator) 

Let D be a connected open set in P d with unit volume. Let P = dx be the Lebesgue 
measure on D. Let T be defined on C™(D) via 

* 3 

where the functions satisfy 

a tl EL 2 oc (D) gLeLUD) i,j,k — 1, . . . ,d, 



then the pre-structure (D, B(D), P, C%°(D), T) is closable. 

Let us now come back to the question of using Dirichlet forms instead of second order 
vectors as germs of semi-martingales. 

The first reason is the closedness property. That gives all the power to this theory. It 
is similar to <T-additivity in probability theory. Without the closedness property, we have 
an apparently more general framework (as additive set functions are more general than 
ex-additive ones), but it becomes impossible to say anything on objects which are defined 
by limits, error propagation is limited to explicit closed formulae. Instead, this closedness 
property allows to extend error calculus to infinite dimensional frameworks and to propa- 
gate errors through typically limit objects as stochastic integrals. As David Hilbert argued 
against intuitionists, more theorems is better. The philosopher Carl Popper made this 
mistake about axiomatization of probability theory emphasing that his system (without 
cr-additivity) was more general than that of Kolmogorov (with cr-additivity). 

What is particularly satisfying is that this closedness property is preserved by prod- 
ucts. Any countable product of error structures is an error structure and the theorem on 
products (cf. [Bou-03b]) gives explicitely the domain of the new T operator. Starting with 
the Ornstein-Uhlenbeck structure in dimension one, the infinite product of this structure 
by itself gives the Ornstein-Uhlenbeck structure on the Wiener space. Less surprisingly, 
the image of an error structure, defined in the most natural way, is still an error structure, 
as an image of a probability space by a measurable map is still a probability space. 

The second reason is related to simplicity. Let us come back to the first remark at the 
beginning of this part. We said that starting with a centered error, centeredness is lost 
after a non linear map. But what is preserved by image? Which property is an invariant ? 
It is the global property of symmetry with respect to a measure. If the operators describ- 
ing the error are symmetric with respect to some measure, the image of the error has still 
this symmetry with respect to the image measure. Centeredness is nothing but symmetry 
with respect to Lebesgue measure (not a probability measure, a a-finite measure but this 
doesn't matter really here). 

The gradient and the sharp (#). 

In addition to the operators T and A we will need the notion of gradient which is a 
linear (Hilbert valued) version of the standard deviation of the error. 
Definition. Let H be a Hilbert space. A linear operator D from ID into L 2 (JP,H) is said 
to be a gradient (for S) if 

\/u G ID r>] =< Du, Du > H . 

A gradient always exists as soon the space ID is separable. It satisfies necessarily the 
chain rule : 

Proposition Let D be a gradient for S with values in H. Then Vw G E) n ,VF G C 1 fl 
Lip(TR n ), 

n dF 

D[F o u] — ^2 — o uD[ui] a.e. 

i=l " X i 

What we denote by the sharp # is a special case of gradient operator when 7i is chosen 
to be L 2 (n,A,F) where (Q,A,TP) is a copy of *A, IP). It is particularly usefull for 



Let us give some definitions and notation we will need later on about the weighted 
Ornstein-Uhlenbeck structure on the Wiener space : let B be a standard Brownian motion 
constructed as coordinates of the space C([0, 1]) equipped with the Wiener measure and 
let a be a positive function in L[ oc [0, 1], there exists an error structure (cf. [Bou-03b]) 
satisfying 

T[/ u{s)dB s ] = f a{s)u 2 {s)ds 
Jo Jo 

for u G C([0, 1]). It is the mathematical expression of the following perturbation of the 
Brownian path : 



u(s)= I 8 dB u ^ ( S e-^ £ dB u + f S ^Jl-e-^dB u , 
Jo Jo Jo 

where B is an independent standard Brownian motion. This structure possesses the 
following ^-operator : 

(£ u(s)dB s )* = jf 1 y/a(^)u(s)dB a , Wu G L 2 ([0, 1], (1 + a)dt), 
which satisfies for regular adapted processes H 

(f 1 H s dB s )*= f 1 J^s)H s dB s + C HfdB s . 
Jo Jo v Jo 

Let us end this part by a comment on the passage from a random walk to the Brownian 
motion in the context of erroneous quantities. Donsker's theorem says that if U n are i.i.d. 
square integrable centered random variables, the linear interpolation of the random walk 
YJk=i Uk i-e. the process 



1 



[at] 



X n (t) = -= £ tf fc + (ra* - [nt])U [nt]+1 



>fc=i 



for t G [0, 1], where [x] denotes the entire part of x, converges in law on the space C([0, 1]) 
equipped by the uniform norm to a Brownian motion. Invariance principles follow giving 
a way to approximate properties of the Brownian motion by the corresponding ones of 
the random walk. A quite natural question is how this may be extended to the case 
where the U n are erroneous. To extend weak convergence of probability measures we 
use convergence of Dirichlet forms on Lipschitz and C 1 functions. Then supposing the 
errors on the C/ n 's are equidistributed and uncorrelated, the error structure of the process 
X n converges to the Ornstein-Uhlenbeck structure on the Wiener space (cf. [Bou-05a]). 
Invariance principles follow giving approximations of the variance of the error of Brownian 
functionals, for example for the sup-norm of the paths : 

Er[||X n (f)||oo] =Er[-L max \S k \] - E[ [' (D s [\\.\U) 2 ds] = E[T] 

\/n l<k<n Jo 



where D denotes the Ornstein-Uhlenbeck gradient with values in L 2 ([0, 1]) and T is the 
random time where the absolute value of the Brownian path reaches its maximum. 

III. Propagation of the error due to the Euler scheme. 



If an asset X is represented by the solution of an sde, prices of options, hedging port- 



discovered by Kurtz and Protter in the spirit of a functional central limit theorem takes 
the form of a process solution to an other sde. In order to propagate this assymptotic 
error through stochastic calculus, we have to take the derivative in a suitable sense of non 
differentiable functionals as stochastic integrals. This may be performed by the theory of 
Dirichlet forms. Let us recall the situation. 

The error due to the Euler scheme. 

In 1991 Thomas Kurtz and Philipp Protter obtained an asymptotic estimate in law 
for the error due to the Euler scheme 

In the simplest case, considering the sde 

X t = x + f a{X s )dB s + f b(X s )ds, 
Jo Jo 

if X™ is the Euler approximation of X t and U n = X n — X then (B, ^JnU n ) converges in 
law to (.£>, U) where U is solution to the linear sde 

dU t = a\X t )U t dB t + b'(X t )U t dt + -j=a'(X t )a(X t )dW t , U = 0, 

where If is a Brownian motion independent of B. 

Such an "extra-Brownian motion" appeared in a work of H. Rootzen [Roo-80] who 
studies limits of integrals of the form J ip n (s)dB s where ip n is an adapted process. In the 
case where J * f(B s , s)dB s is computed by the Euler scheme 

/ ij n {s)dB s = Y,f{B^ ) i/n){B 1±1 -B^) + f{B m) [nt}/n){B t -B^) 

JO ~q n n n n n 

he obtains for regular / 

Qf i> n dB - J q f(B„ s)dB s ^j A ±= J o f' x (B s , s)dW s . 

This kind of result is restricted to adapted approximations. As Wong and Zakai have 
shown (1965 ) other natural approximations of the brownian motion give rise to stochastic 
integrals in the sense of Stratonowitch 

The discovery of the asymptotic error due to the Euler scheme has been followed 
by a series of works which extend it to the case of an sde with respect to a continuous 
or discontinuous semi-martingale and which obtain some statements as necessary and 
sufficient conditions ([Ja-Pr-98], [Ja-Ja-Me-03]). 

In addition, asymptotic expansions have been recently obtained by the stochastic 
calculus of variation [Ma-Th-03] . 

In the sequel, we shall consider the result of Kurtz-Protter in dimension 1 under the 
following form: 

Let X t be the solution starting at x to the sde 

dX t = a(X t , t)dB t + b(X t , t)dt, 
let X™ be the approximate solution obtained by the Euler method, which may be written 

X " = x ; dX? = a(XU, [nt]/n)dB t + b(XU, [nt]/n)dt 

n n 

and let = X™ — X t be the approximation error, then if a and b are C 1 with linear 



where the process U may be represented as 

U = dU t = a' x (X t , t)U t dB t + b' x (X t , t)U t dt + ^a' x (X t , t)a(X t , t)dW t 

which is solved by the usual method of variation of the constant : introducing the process 

ft i rt r t 



M t = exp | jT <(X S , s)cLB s - - jf < 2 (X S , s)ds + jf b' x (X s , s)dt 



gives 



Jo 



rt a{X s ,s)a' x (X s , s) 
V2M S 

Let us consider the weighted Ornstein-Uhlenbeck error structure on the Wiener space 
with weight a as explain above. If the coefficients a and b are regular, then X t G ID and 
Xf satisfies 



(*) X* = fa' x {X s , s)XfdB s + fa(X s , s)J^(s)dB s + l\ x {X s , s)Xfds 
Jo Jo v Jo 

Comparing with the equation of the asymptotic error due to the Euler scheme 

(**) U t = fa' x {X s , s)U s dB s + fa(X s , S ) °^' 8 W , + fb' x (X s , s)U s ds 
Jo Jo v2 Jo 



shows that 

- if we could take a random and adapted weig ht a(t) = \a'^{X u t) 

- if the obtained structure is closable with carre du champ and if the calculus of the 
^-operator is still (*) 

then X* would be the asymptotic error due to the Euler scheme, and we would be able 
to propagate this error through the stochastic computations obtaining the variance of the 
error on any r. v. Y E ID by the equation T[Y] = 1E[F# 2 ]. 

The Ornstein-Uhlenbeck structure with random weight. 

From now on a is a measurable random process defined on the Wiener space, non 
negative, non necessarily adapted. We assume that this process satisfies IE Jq 1 a t dt < +oo, 
and a(iv,t) > k(t) > IP x dt-a.e. where k is deterministic. 

Let us denote ID^ U the domain of the Ornstein-Uhlenbeck structure with deterministic 
weight k and D k ou its gradient. On the domain 



ID = \y e JD k ou : £ n(D k AYKt)f^}dt < +oo 
which is dense, the form 



kit) 

is Dirichlet and admits 



as carre du champ operator. 

Indeed, let V be the space of linear combinations of exponentials of the form Y 



Let X n be a Cauchy sequence in L 2 and for £. Let X be the limit of X n in L 2 . Then 
X n is Cauchy for £ k u which is closed, hence X G ID^ and there exists a sub-sequence AV 
such that 

D k a [X n ,] D k ou [X] Exdt-p.s. 

and by Fatou's lemma 

= jf' E[lim(B;„[X,„]) 2 ^|]<i« < liminf jf' E[(Z£jX„.]) 2 |j|]<ft < +«d 

since X„ is Cauchy for £. Hence X G ID. Now again by the Fatou's lemma we show as 
classically that X n converges to X in D. 

Contractions operate on (£, ID) by the functional calculus for D k u hence (£, E>) is a 
Dirichlet form. The definition of the carre du champ operator (def 4.1.2 of [Bo-Hi-91]) is 
satisfied. 

The generator (A, VA) is given by 

VA= {FgE 3GgL 2 MH G ID ^E^ 1 £&[F]l£,[tf]£gdt = - < G, >} 
AF = G 

hence if F G £>A then f[|r/jF] G dom<^ and 

where 5 k u is the Skorokhod integral with weight k. 
Adapted case. 

Let us now add the hypothesis that a is adapted. If h is in L°°(IR + ) 

Er[F, f 1 hdB]] = E[F f 1 h(s)a(s)dB s }. 
Jo Jo 

If F, G G ID n L°° 

E[G < DF, >] = — E[F < DG, h^a>] + E[FG ^ hadB}. 
And if f is adapted and in dom5 

5[v]= [ v sy /c^dB s . 
Jo 

At last, for finance, the following properties are important, they use the fact that a is 
adapted 

ApE[X|^J] =nA s [X]\F s ] 
where A s is constructed as A with the weight a(t)l{ t < s }, 

D\E[X\F,]](t) = 1E[D[X}(t)l t < s \F s } 

E[. \F S ] is an orthogonal projector in ID 



where # s is constructed as # with the weight a(t)l{ t < s }. If X is jF r measurable, then 
AX, T[X] are jF r measurables. 

Concerning the operator # we have the formulae 

Qf 1 e s * = £ c # ^ s + £ CsV^ s dB s 

Hence formula (*) is satisfied. 

Application to diffusion models. 

Let us consider the following model of an asset 

dX t = X t a(X t , t)dB t + X t r(t)dt 
and let us put on the Wiener space the Ornstein-Uhlenbeck structure with weight 

_ a' 2 (X t ,t) _ (a(X t ,t) + X t a' x {X u t)f 
at ~ ~~2 " 2 

which represents the asymptotic error due to the Euler scheme, a is supposed to be 
strictly positive, C x and Lipschitz and the preceding hypotheses on a are assumed. 

Such a modelisation is coherent. The error is attached to the asset X and any func- 
tional of X, including the Brownian motion itself and its error may be computed thanks 
to the equation 

which gives 

(B t )* = f J^(sjdB s T[B t ]= f a(s)ds. 

Jo v Jo 

Let us show how financial calculi may be performed before proposing some comments on 
the use of such an analysis. Puting M t = exp{ J * ^fa~ s dB s — \Jq a s ds + J l r (s)ds} we have 

T[X S , X t ] = M s Mt I u \y 2 U ' } a u du. 
Jo M£ 

The price of a European option with payoff f(X T ) at exercise time T 

V t = n(eW-f t T r(s)ds)f(X T )\F t ] 



becomes erroneous (in the sense of error structures) with an error obtained thanks to the 
#: 



T 



T[V t ] = (exp-2/ r{s)ds){nf{X T )M T \F t }f T [X,] 



M 2 

T[V s ,V t ] = (eM-Jj{u)du-£r{v)dv)^ E[/' {X T )M T \F s ]E,[f (X T )M T \^ 
The quantity of asset in the hedging portfolio is 



1 



and we have 



T[H t ] = (e W -2j t r(s)ds)(n^(f(X T )M T + f(X T )Zj)\T t }) 



M? 



with 



Zj = [ L s dB s - [ ^cT s L s M s ds 



Jt Jt 



L s = a'* 2 (X s , s) = 2a' x (X s , s) + X s a^(X s , s) 



It is still true, as in the case of deterministic weight (cf [Bou-03b]), that the proportional 
error on X t divided by the volatility : 



is a finite variation process (cf [Ba-Ma-Ma-Re-Th-03]) on the "feed back" effect). 
Discussion. 

Thanks to this construction of an error structure, i.e. a local Dirichlet form with 
squared field operator, on the Wiener space, hence by image, on C([0, 1]) equipped with 
the law of the process X, we have at our disposal a powerful mean to propagate the 
error done on X toward sufficiently smooth functionals of X. In order to assess the 
interest of this tool, the question arises of knowing whether the propagated error is the 
same as the one we would obtain by a direct computation of the functional thanks to the 
approximation X n of X. For instance, in the simplest case, does the convergence in law 



hold for / G C 1 fl Lip ? Can we justify an asymptotic calculus principle which says that 
the Dirichlet form allows effectively to compute the errors on the quantities which are 
erroneous because of the approximation X n of X ? We will not exhaustively examine 
this principle here, for it is a too large enterprise. Nevertheless, in the important current 
of research whose fruitfulness has been confirmed these last twenty years, which may be 
called the "tightness programm" , the authors, among which we must at least quote P.-A. 
Meyer, W. A. Zheng, J. Jacod, A. N. Shiryaev, A. Jakubowski, J. Memin, G. Pages, 
T. G. Kurtz, P. Protter, L. Slomihski, D. Talay, V. Bally, A. Kohatsu-Higa and many 
others, have already done a major part of the work by stating their results of convergence 
in law, of stable convergence, of tightness of processes, under a sufficiently general form 
for propagating iteratively the properties through stochastic integrals and sde's in the 
semi-martingale framework. 

Let us give some results in the direction of this asymptotic calculus principle keeping 
the hypotheses of the present part III. 

Let F be a real function of class C l and Lipschitz defined on C([0, 1]) equipped with 
the uniform norm. Such a function satisfies 



where the mapping x \— > F'(x) is continuous and bounded with values in the Banach space 
of Radon measures on [0, 1], e x {h) is bounded in x and h, and goes to zero when h — > 
in C([0, 1]). Then we have 



X t ' a{X t ,t) 



F(x + h) = F(x)+ < F'(x), h > +\\h\\e x (h) Vx, h e C([0, 1]) 



Proof. The equality in the right hand side comes from the functional calculus in error 
structures (see [Bou-05a]). Puting U n = X n - X as before, the fact that ^\\U n \\e x {U n ) 
tends to zero in probability, reduces the proof to the study of the convergence in law of 

< F'(X), V^U n >=V^j W - X t ) F\X){dt) 

to fxf F'(X)(dt). Considering the measure F'(X)(dt) as the differential of a finite 
variation process adapted to the constant filtration Q t = £>(C(]0, 1])), the fact that the 
process to be integrated y/nU n converges stably to X* implies (cf [Ku-Pr-91b] thm 2.2) 
that the stochastic integral / ^JnU n F'(X)(dt) converges in law to fXf F'(X)(dt). 

We obtain also the convergence in law of the stochastic integrals H.^/nU n 4- H.X* 
for H deterministic or adapted and that of 

M jf 1 / w, S )dx: - jf 1 f(x s , S )dx s ) 

to 

( /' f(X s , s)dX s )* = f 1 f'(X„ s)XfdX s + f 1 f(X s , s)dXf 
Jo Jo Jo 

for / C l and Lipschitz. 

More generally, we can make more explicit the research programm of determining the 
domain of the asymptotic calculus. 

Let X n and X be two random variables with values in a measurable set (E,!F), and 
let a n be a sequence of positive numbers. Let V denote a set of simple functions included 
in L 2 (W X ) and in L 2 (W Xn ) Vra. Let us suppose that there exists an error structure 

s = (E,r,w x ,-D,r) 

such that V C ID and Vip G T> 

(1) hma n lE[(v9(X n ) - V {X)f] = E[r [</?]] 

we shall say that the asymptotic calculus principle extends to V for V C V C ID if the 
limit (1) extends to i[> G T>. 

If, as above, a ^-operator is available (which occurs as soon as ID is separable), in 
order to prove (1) on D, since # is a closed operator, it suffices for any xjj G T> to find a 
sequence ip p G T> such that 

i) ^ p ^^inL 2 (P x ) _ 
m) <y9* converges in L 2 (JP X x Px) 

iii) a n !E[^(X n ) — -^(X)) 2 ] may be approximated uniformly in n by 

a n nMXn) - MX)) 2 )- 

When (E, T) is a normed vectorspace, obtaining (1) from a convergence in law of y^a^((p(X n ) — 
<p(X)) uses generally a uniform integrability of a„||X n — X\\ 2 . We shall go deaper in this 
problem in a separate work. 

Let us end by some remarks from the point of view of finance. The interest of consider- 
ing a financial asset as erroneous is not evident since it is one of the best known quantities 

„ — i-; 4- — 1 „ o,,„u u„ f„c TT3„,, noul\ 



bid-ask. But this would rather justify specifically constructed error structures instead of 
the one induced by the Euler scheme. This error structure is relevant only in order to 
assess the errors in Monte Carlo simulations performed to calculate financial quantities 
in a given model. 

Several authors ([Du-Pr-89], [Ha-My-05]) remarked that the stochastic integral which 
is the active hedge of a future contingent claim, in a model where the underlying asset is 
a semi-martingale, is an instance of application of limit theorems on discretization errors. 
This is different from the Euler scheme error and it would be worth to examine this error 
from the point of view of an asymptotic Dirichlet form. 

A more general and complete study of the bias operators and the Dirichlet form yielded 
by an approximation, with applications related to the part III of the present paper, is to 
appear, ([Bou-06]). 
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